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Abstract 

We present a multiscale description of hydrodynamic turbulence in incompressible 
fluid based on a continuous wavelet transform (CWT) and a stochastic hydrodynamics 
formalism. Defining the stirring random force by the correlation function of its wavelet 
components, we achieve the cancellation of loop divergences in the stochastic pertur- 
bation expansion. An extra contribution to the energy transfer from large to smaller 
scales is considered. It is shown that the Kolmogorov hypotheses are naturally refor- 
mulated in multiscale formalism. The multiscale perturbation theory and statistical 
closures based on the wavelet decomposition are constructed. 

PACS: 47.27.Gs, U.lO.Gh 

1 Introduction 

The statistical description of a fully developed hydrodynamic turbulence is based on the 
Kolmogorov hypotheses PP on the self-similarity of the velocity fluctuations of different scales. 
However, the Kolmogorov-Obukhov analysis |21 IHj does not provide a rigorous mathematical 
definition of the "fluctuation of scale /". In the literature this is tacitly understood as 
Fourier components with wavenumbers approximately equal to the inverse scale ^ ~ 7^ 
the analysis is performed in wavenumber space. This definition meets global characteristics 
of the fully developed isotropic turbulence, but, being based on the Fourier transform, is 
essentially nonlocal and therefore hardly applicable to such important properties of the fully 
developed turbulence as the coherent structure formation. To catch the local properties of 
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the turbulent velocity field, the decomposition into localized wave packets and wavelets have 
been performed by many authors [U El El 13 El E] • 

Most of the wavelet applications to turbulence are restricted either to an analysis of the 
measured turbulent fields with the "wavelet microscope", capable of simultaneous analysis 
of the same velocity field at different resolution [71 El El ? or to a numerical solution of the 
Navier-Stokes equations (NSE) in the wavelet basis [Ul jTUl HH 112] • The application of the 
continuous wavelet transform with the derivatives of the Gaussian taken as basic wavelets 
was already used in the analytical study of the NSE ^3 , in particular, in research on 
energy dissipation. However, at least to the authors knowledge the wavelet decomposition 
has not been yet applied to the stochastic iterative solutions of the Navier-Stokes equation or 
in the framework of the field theory approach to statistical hydrodynamics ITTj. 

Our interest in extending the stochastic hydrodynamics approach by wavelet-defined 
random processes is stimulated by recent developments in application of the field theory 
and renormalization group methods to the fully developed isotropic turbulence, see e.g. 
[TH] for a short summary. The "field theoretic methods" are taken to mean the stochastic 
diagram technique, the functional integral representation of the characteristic functional and 
the renormalization group methods. The latter, being inherited from the theory of critical 
behavior in equilibrium phase transitions clearly demonstrates the need for a scale-dependent 
random measure in the field theoretic approach to turbulence. Such measures, having been 
already in use for phenomenological studies of multifractal behavior and itermittency ^Hl HI 
E], have not been yet used in analytical approach of stochastic hydrodynamics. 

The aim of this paper is to extend the wavelet representation of the stochastic Navier- 
Stokes equations in such a way that the probability distribution of the stirring force, used to 
compensate the energy dissipation, is defined for the wavelet coefficients of random force (i.e. 
for the scale components of forcing). In our description the velocity field wavelet coefficients 
ui{x) attain the Kolmogorov meaning of local velocity fiuctuations of scale / at a given point 
X (A particular case of the difference of two Gaussians or delta-functions was considered in 
the literature 01201 )• Defining the random force in the space of wavelet coefficients we get 
an extra analytical fiexibility: there are random processes in the space of wavelet coefficients 
with different correlation functions whose images under inverse WT coincide in the space of 
common random functions. Tuning the random force correlation function in the space of 
wavelet coefficients, in stochastic hydrodynamics formalism, we get rid of loop divergences; 
for a special type of narrow-band forcing the contributions to the response and the correlation 
functions are explicitly calculated. 

It is shown that the Kolmogorov hypotheses, the statistical closures of moment equations, 
the stochastic hydrodynamics approach and the Wyld diagram technique are naturally refor- 
mulated in multiscale (wavelet) formalism. Besides, the consideration of random processes 
depending on scale explicitly {ui{x)ui'{x')) = C{x,x',l,l') gives a possibility of the pertur- 
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bation expansion converging without introducing an ultra-violet (UV) cutoff wavenumber, 
and an extra contribution to the Kolmogorov energy dissipation term {uf / 1) . 

The rest of this paper is organized as follows. In Section 2 we review the methods of 
solution of the NSE by the Fourier decomposition and by the wavelet decomposition. In 
Section 3 the stochastic hydrodynamics approach is reformulated for the random processes 
explicitly depending on scale. A regularization of the perturbation expansion for a random 
force acting at a single scale is presented. Energy dissipation rate and energy flux in the 
multiscale formalism are considered in Section 4- Section 5 presents a generalized form of 
the Kolmogorov hypotheses formulated in a multiscale framework. In Conclusion we discuss 
some perspectives of the method. 

2 Hydrodynamics of incompressible fluid in wavelet 
representation 

In analytical studies of the hydrodynamic turbulence the basic role is played by the Navier- 
Stokes equations: 



Our consideration of hydrodynamic turbulence, based on the Navier-Stokes equations, as- 
sumes a fully developed homogeneous isotropic turbulence far from any boundaries. Al- 
though the adequacy of the stochastic NSE to fully developed turbulence still remains an 
open problem, a significant progress has been achieved in studying simple models, proving 
the Kolmogorov spectrum from basic principles |2I] and calculating the anomalous scaling 
corrections [22 1231 121] • Turbulence in an incompressible fluid considered in these settings is 
more frequently studied in a wavenumber space rather than in a real space. The advantage 
of wavenumber space is a very simple form of the Laplacian - this enables to eliminate the 
pressure term from the NSE. The price paid for this simplification is the non-locality of the 
Fourier transform that completely hides all information related to the real space distribu- 
tion of the velocity field. The discrete Fourier transform used in numerical simulations also 
imposes periodicity on the system. The usefulness of the wavenumber representation stems 
from two basic reasons: the existing of fast FFT algorithms and pseudospectral methods 
[23 I2n] and the direct experimental interpretation of the power spectra density of velocity 
fluctuations [2ZI. For the incompressible fluid the pressure term can be eliminated from the 
NSE ([!} by the substitution 






(2) 
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In the wavenumber space representation 

«.(x,t) = J e^'^-«.(k,t)^, u,{^,t) = J e<")t..(k,a;)^|;^, (3) 

both the inverse Laplacian operator in (j21) and the incompressibihty condition V-u = 
are simphfied and the NSE system becomes a system of integro-differential equations 

{dt + iyk^)uiik,t) 

where Mijk(k) 

The system of equations (jH) allows for perturbative calculations, statistical closures of 
moment equations etc., but is not local in the coordinate space and therefore is of little use 
when studying the effects locally produced by fluctuations. The system of equations (jD) is 
complete and, being correctly solved numerically, gives reliable results. This however requires 
taking into account a tremendous number of Fourier modes, that is hardly bearable even 
for modern supercomputers; alternatively certain block-averaging procedures in wavenumber 
space can be applied, e.g. the spectral reduction j26j. 

The idea of studying turbulence using the space-scale or wavenumber- scale representation 
is not a new one. The band-pass filtering of the turbulent signals was used to study inter- 
mittency long ago [2H]- Later the prototype of wavelet cascade model was given in The 
wavelet transform (WT) is known to be an excellent local tool, widely used in a turbulence 
data analysis ^ and numerical simulations IHIHllIll- WT in real space, as a particular type 
of filtering j^H] , reveals the formation of the coherent structures jH] and is useful to study the 
local energy dissipation effects related to filament formation ^1]. WT in (a, k) space enables 
to enhance the idea of band-pass filtering by separating the contributions of different scales. 

In this paper we use continuous wavelet transform (CWT) to derive the equations for the 
fluctuations of different scales. This analytical consideration goes along with the filtering ap- 
proach fl^ IHU] , when considering moment closures for a fully developed isotropic turbulence. 
The crux of our approach is the application of continuous wavelet transform in the stochastic 
hydrodynamics framework: the WT is applied to both the velocity field and the random force 
stirring the turbulence. Working in the space of wavelet coefficients, rather than in a real 
space of velocity fields, we can define the random stirring force, which is essentially used in 
RG calculations, see e.g. EI - in a way that allows for getting rid of loop divergences 
in the stochastic perturbation expansion of the velocity field statistical momenta |32j . 

For simplicity we restrict ourselves with the homogeneous isotropic turbulence and the 
isotropic wavelets iplx.) = In this case, the wavelet transform of the velocity field 



d q 



^Mijk(k)uj{ci, t)uk{'k - q, t). 



(4) 



-- [kjPik{k) + kkPij{k)] , Pij (k) = 5ij - 
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u{x,t), taken with respect to the basic wavelet and the corresponding reconstruction 

formula are 



u„(b,t) = / A^^(^-^)u(x,t)d'^x, (5) 



a 



X 



b\ , dad^h 



We perform the wavelet transform only in the spatial argument of the velocity field because 
we need the spatial resolution. Using the norm instead of L^, we provide wavelet coeffi- 
cients Ua{h,t) with the same {LT~^) dimension as the velocity field u{x) itself. The wavelet 
coefficients Ua(b,t) are referred to hereafter, as the components of the velocity field corre- 
sponding to scale a; ip{x) is referred to as an analyzing function used to measure the scale 
components. 

For practical calculations it is often convenient to express wavelet transform ()5pfi|) in (a, k) 
representation taking the Fourier transform in the spatial argument. In Fourier form the 
direct and inverse WT are: 

Ua(/i;) = ^(ak)u(A;), u(A;) = /" ^?/^(ak)ua(A;), (7) 
where u{k) = u(k, u) is the Fourier transform of the velocity field 



(2vr) 



d+l ■ 



(The Minkovski-like notation x = {:K,t),k = (k, cu) is used.) Therefore, the wavelet trans- 
form (0) can be considered as a frequency filter that conveys the harmonics with typical 
wavenumbers of order ^ and is localized close to point b. 

The only restriction imposed on the basic wavelet ip to make the wavelet transform 
invertible - the admissibility condition - is the finiteness of the normalization constant C^: 

= / ' M^rfa < oo. (8) 



For a real- valued basic wavelet V'(x) we can restrict the integration to the positive frequencies 
only 
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If the basic wavelet is also isotropic ipi^x.) = '?/'(|x|), we get 

where 5*^ is the area of the unit sphere in R'^. 

In this paper we will assume a G M+ integration and the isotropic real wavelets. Thus, 
the decomposition of the velocity field with respect to the basic wavelet ip takes the form 

«(x,t) = T^r^Tl/ d'^h^(^)n^{h,t) (10) 
Cv. Jo jRd \ a J 

.(x,t) = ^r^/ |^e<---)^(akK(k,c.). (11) 

We drop the integration limits ^ hereafter. 

Substituting the wavelet transform (fTTj) into the system of the component equations 
we yield the system of equations for the scale components u.ai{k): 



{q)Ua2k{k-q) 



dai da2 d'^^^q 
ai a2 (27r)'^+i 



M?;",^"^(k,q,k-q) = 7A(ak)M,,fc(k)^(aiq)V'(a2(k-q)). (12) 



Let us derive statistical closures for the scale components. For this purpose we take Eq. ^2 
and its complex conjugate 



{^t + uk')Ua^{^^,t) = (^) /M---(k,q,k-qX,(q,tX,(k-q,t)^^0 

[d, + uk^)Ijk;f) = (^)VM^r^(k,q,k-q)<,(q,t')^a.fc(k-q,t')^^0, 

(13) 

multiply the first equation by Uai(k, f), sum up over the vector index i and take the statistical 
averaging ( ). Doing so, we get 

(5, + .k^)j;(UM0^a.(k,t)) = (^)7t^^(^^""^^'^'^-^) 

{Uai{k, t')Uaij{c[, t)Ua^kik-q, t)) . 

Applying the same procedure to the second of the equations (|T3| and summing up the results 
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at coinciding time arguments t = t', we get the moment equation 

{itaiik, t)uaj^j{q, t)ua2k(M-q, t)) + h.c, (14) 

which is different from its plane wave counterpart only by extra scale indexes and extra 
integrations in scale logarithms / — (octaves). To express the third order moments in (fT^ 
via the second moments, we must substitute 



Un 



(k,t) = (^) /lrf.Gr(k,t-.)M---(k,q,k-q) ^^^^ 

'"aij(q, ■S)Ma2fc(k-q, S)^^^j^, 

where G";'^''(k, t — s) is the response function. The difference from the standard plane- wave 
approach |H] is that additionally to the summation over vector indices we have to sum up 
over octaves to integrate over J ^ in each scale variable. So the statistical closures can be 
reproduced for the scale components. 

In the zero-th order approximation (with no interaction term: M(-) 0) the bare 
response function is given by 

G;r^^"(k,t-.) = 6, r Si^-^oH ,~^^it~s)p ^ ^^^^(^ _ a,)aoe-^^'\^-^^. (16) 

The full response function, in view of the component equations fll3|) . satisfies the integro- 
differential equation 



(9,+z/k2)Gr (k,t-.) = 2 J M---(k,q,k-<i)u.,,(q,t)Grnk-<i,t-.) 



The substitution of (fT^ into (fT^ gives a relation between the second and the forth order 
moments of the scale components: 



{dt + 2uk') J2 (^-(k, t)uai{k, t)) = 2 (^^^ j 



dai da2 d'^ki da^ da^ d'^k2 



^M;}'^,i'^^(k,ki,k-ki)(MM^ia,,(ki,t) f rfsGr°(k-ki,t,s) 

«0 J-00 

Mf;;^'^^(k-ki, k2, k - ki - k2)n,3,(k2, sX/(k-ki-k2, s)) + h.c. (17) 

The forth order moments {uuuu) can be further decomposed into the sum of all pairs 
{uu) {uu) using a stochastic perturbation expansion. 
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3 Stochastic hydrodynamics with multiscale forcing 



The stochastic hydrodynamics approach consists in introducing random force in the Navier- 
Stokes equations and calculating the velocity field momenta (u(xi) . . . u(x„)) using the 
stochastic perturbation theory, pioneered by Wyld or the functional integral formal- 
ism. 

In a coordinate representation, the stochastic NSE is written in the form 

+ {u-V)u = uAu-V +r]{^,t). (18) 



dt ' ' \p 

The random force correlator {r]i{x)rij{x')) = Dij{x — x') should obey certain conditions to 
make the resulting theory physically feasible and the perturbation expansion suitable for 
analytical evaluation. First, the energy injection by random force should be equal to the 
energy dissipation; secondly, the forcing should be essentially infrared (IR), i.e. it should 
be localized at large scales; third, it is desirable to have a parameter to control the IR 
divergences (when the size of the system tends to infinity). 

To exclude the pressure term in (fT^ using the incompressibility condition, the Fourier 
representation is used 

[-IUJ + z/k^) Mi(k, ^) - J ■^^:^j:p^^iifc(k)Mj(g)Mfc(/^ -q) = Viik). (19) 

The random force correlator is usually taken in the form 

(UhMh)) = {2TxY+H''+\k, + A;2)P,,(ki)/^(ki), (20) 

where the function -D(|k|) has a suitable power-law behavior. In the simplest, but not very 
feasible physically, case D{\s) = Dq = const, we deal with the white noise, (5-correlated in 
both space and time. 

What is most realistic physically, is to have a random force concentrated in a limited 
domain in fc-space Amj„< |k| <Amax- This case, however, is difficult to evaluate analytically 
in the perturbation theory jT^I. In the multiscale approach we are going to present, we solve 
this problem by constructing a noise acting in a limited domain of scales a in (a, k) space. 

In the (a, k) representation, having excluded the pressure by standard means of the 
orthogonal projector, the Eq. ()18|1 leads to a system of integro-differential equations for the 
scale components Uai{k): 

{-lU + ■u\i^)Uai{k) = flai{k) (21) 



+ y"il^'^;'^^(k,q,k-qX,(gX,(A;-g). 



ai a2 (27r)^+i 
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Now we face the problem of appropriate choice of the force correlator 

It was shown in the previous paper [32J that the (a, k) representation provides an extra 
analytical flexibility in constructing random processes with desired correlation properties in 
a coordinate space. For instance, the random process given by wavelet coefficients with the 
correlation function 

{fjaMfjaM) = (27r)'^^5'^(ki + k2)ai5(ai - 02)^0 (22) 

possesses the same correlation properties in Mf^ coordinate space as the white noise does. 
Casting r]{k) in terms of fja{k) by means of (|7j) we get 

J J ai a2 

= {27cy^6'{k, + k2) [ ^^V^(aiki)^(a2k2)ai5(ai-a2) 

J Oi 02 

= (27r)'^5'^(ki + k2)Z^o 

that coincides with the correlation function of the white noise. The direct wavelet transform 
of the white noise ri{x.) ^fi(k) ^fia{k) apparently leads to another result 

(fiaMfjaM) = (27r)'^5"(ki + k2)Z^oV^(aiki)^(a2k2), (23) 

which is different from (j^^ and explicitly depends on the basic wavelet ip. 

Physically, the scale-dependent processes obeying and pHjl . respectively, describe 
quite different processes: fluctuations of the former type ()22|) are mutually correlated only 
for coinciding scales (01 = 02), while for the latter case all fluctuations are correlated. 

Exactly as in the standard wavenumber space approach fT^ I16[ IT7j. we can generalize 
the 5-correlated force (j22j) assuming its variance to be dependent on both the scale the 
wave vector: Dq — > D{a,k.). Taking into account we deal with the incompressible fluid in d 
dimensions, we can put down a general form of the desired force correlator 

{VaAki)Va,Ak2)) = {2T,f+^5''+\k, + k2)^a^5{a, - a2)P.,(kO/^(ai, |ki|). (24) 

The (5-correlated random force in the wavenumber space does not provide an adequate de- 
scription of hydrodynamic turbulence, for it gives an energy injection in all scales, small 
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and large. In physical settings, the fluid is usually stirred at a predetermined scale, or in a 
narrow range of scales, comparable to the size of the system. As a simplest model of such a 
forcing, we can consider a force acting on a single scale by choosing 

D{a,k) = Doao6{a — ao). (25) 

Now let us turn to the perturbative calculations. The stochastic diagram techniques for 
the component flelds Ua{k) stems from Eq. (j2ip and is a straightforward generalization of 
the Wyld diagram technique for the Fourier components u(k): 

ua^ik) = Go{k)Uk) + Go{k)(^y f^^-f^ (26) 

M^;i'^^{k, q, k - q)UauiQ)i^aAk - q), 

where Go{k) = {—luj + z/k^)^^ is the bare response function for Fourier component. To keep 
the scales and wavevectors on the same footing and make the notation covariant in that 
sense, we can rewrite using the response functions bearing scale indices explicitly (fT^. 
Thus, we get 

. / . (A) Y p.) 

M,7,^"^(k, q, k - Ci)ua,M)i^aAk - q) 

The Feynman expansion for the scale component flelds Uai{k) can be derived either from 
or (j2Zj). Iterating the Eq. once, we get the one-loop contribution to the response 
function: 

Uai{k) = Go{k)7]a^{k)+Go{k)[—] — — jTT^ (28) 



G^J J ai a2 (27r)^+i 
M;5'^,i'^^(k,ki,k-kiX,(fci)[Go(A;-A;i)w-(A;-A;i) 

Ua3l{k2)Ua^m{k- ki-k2) 

As usual, we assume the random force to be gaussian, with all odd correlators {rji . . . ?72fc+i) 
vanish identically. 



10 



Following ^H] ; we introduce a formal parameter of the perturbation expansion A in the 
interaction term (M^^j^^"^ XM^f,}"^). In the final end the initial value A = 1 should be 
restored. The validity of considering A as a small parameter of perturbation expansion is 
justified by renormalization group methods ^3 ^21 ^] ; see also [21] for recent developments 
and generalizations. 

In the zero-th order of perturbation expansion the response function does not depend on 
scale and coincides with Go for all scale components: 

Uaiik) = Go{k)fja^{k). (29) 

In the O(A^) and the next orders of perturbation expansion the standard stochastic diagram 
techniques, used by many authors [23 UHl El, is reproduced with the difference that: (i) 
each vertex, each response and correlation function attain scale superscripts; (ii) integration 
over octaves J — is performed over all pair-matching scale indices. Mathematically, this 
means that each diagram line corresponding to the plane wave component Ui{k) in standard 

techniques now attains an extra wavelet factor and becomes Uai{k) = tp{ak.)ui{k). The 
Feynman graphs and their topological factors of course remain the same. 

Here, for bookkeeping reasons, we present only one loop IPI contributions to the response, 
see Fig. ^ and correlation, see Fig. |2l functions of stochastic hydrodynamics. In the first 
order in the force correlator (7777) (one-loop contribution) we substitute the scale components 
u in the r.h.s. of (|28|) by the zero-th order solutions (j29p . perform necessary scale integrations 
in ai6{a — ct/)^ and average over the random force 

„ „ n\ , ^ n\A\2 I " \ I dai da2 da^ dtti d'^^^ki d'^^^k 



C^p J J CLi (I2 Q-s 0,4 {2'kY+^ (27r)'^+i 
M;^";'^^(k,ki,k-ki)Go(A;i)(r),,,(A;i)r),3i(A;2))Go(A;2)Go(A;-A;i) 
MlZ'^'i^ - ki, k2, k - ki - k2)Go(A;-A;i - A;2)r/,,„(A; -k^- k^). 

The corresponding diagram is shown in Fig. ^ The factor 4 accounts for two possible ways 
to expand the nonlinear term multiplied by two possible ways of taking the random force 
averaging. After substituting the random force correlator ()24j) , the using of the explicit form 
of the interaction vertexes (fT^ . and integrating over the scales, this leads to 

uUk) = Goik)f,Uk) + Ha^)Glik)AX'i—] j,--^ (30) 



C^J J {27rY+^ a, 
|Go(A;i)pA,Kki)Go(A;-A;i)Mfci„(k - k,)iP{a^k)ria,m{k), 
2 \ f dai I 7 ^ , 2 



A,Kki) = P,i{k^){—j J -^ma,k,)\'D{aM- (31) 
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Figure 1: One loop contribution to the response function 
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Figure 2: One loop contribution to the correlation function 



Eq. (jHlj) means that for our special type of the scale-dependent random forcing ()24j] all inter- 
nal parts of the diagrams, that do not carry the scale indices, can be evaluated by substitution 
of the effective force correlator (jHTj) into the standard diagrams drawn in wavenumber space. 
In this way we can easily evaluate the perturbative corrections to the usual response function 
G{k) for Fourier components (see Appendix), and hence evaluate the turbulent corrections 
to viscosity, produced by scale-dependent force. 

Similarly, we can evaluate the contributions to the correlation functions {Uaii{ki)ua2j{k2)) ■ 
In the one-loop approximation, using ()28p and the zero-th order approximation ()29|) . we get 



After the integrations over 04, ag, k^, identically to the response functions calculations, with 
the random force ()24|). we get the one-loop contribution to the correlation function, shown 
in Fig. El 



Ua2j 



= Go{ki)Go{k2){fiaAki)fia2j{k2)) 




M;;^"3'^^(kl,k3,kl-k3)Go(A;3)(r)a3Kfc3)^a4n^(M)G'o(A;4)Go(A;l-A:3) 

{flaAh-ks)fla,nik2-ki))Goik2-k^)Mf^^'''{k2,k,,k2-^A). 




das da^ d'^^^k^ 



a3 as {2txY+^ 



(32) 
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Expanding the wavelet factors ip{a]i) in each vertex and integrating over all matching scale 
arguments, we get the one-loop correction to the correlation function 

|G'o(fc3)|%™(k3)|Go(A;i-A;3)|'Afc„(ki-k3)M,„„(-ki). (33) 

The evaluation of the one-loop contribution to the correlation function is easily performed 
for the above mentioned narrow-band force correlators. The contribution takes the form 

C2(ai, ^2, h) = S'^^^ih + /c2)V^(aiki)V^(-a2k2)Ce//(fci) 

C^ff{k) = 2X'\Go{k)\' J ^^j^A(q)A(k-q)c2(k,q)|Go(g)nGo(fc-g)r, (34) 

where the trace of the one- loop tensor structure C2(k, q) is given in Appendix. The integration 
in frequency argument in the limit of zero frequency (/cq— >0), is not different from that in 
stochastic hydrodynamics in the wavenumber space and gives 



Cffik) = \^\Go{0,k)\'Sd-i J 
(1 -cos2 0)P 



q^-^dqde sin^-' 0A(q) A(k - q) 
z/3g2 ( A;2 -2kqcos9 + q^){k'^ -2kqcos9 + 2q^) 

-[k'^{d-l)-2kqdcos9 + 2q'^{d + 2cos'^9-2)]. (35) 



4(A;2 - 2kq cos9 + q' 

As an example, let us present the one-loop contribution to the effective pair correlator (j34p . 
calculated for the case of single-scale forcing ()25|) with the basic wavelets from the gaussian 
vanishing momenta wavelet family 

^„(fc) = (27r)f(-zk)"exp(-kV2), = (27r)'^r(n), (36) 

with r(x) being the Eulerian gamma-function. 

For the single-scale forcing ()25|) this gives an effective force correlator in the wavenumber 
space 

A„(g) = -^(aog)2"e-('^«'')^ (37) 
r{n) 

Straightforward calculation leads to: 

2A;2|G'o(A;)p5d_ia^"D2 f g^n^p _ 2fcg cos + g2)«e-"o('='-2fc9cose+2g2) 



"^^^ ' 4z/3r(n)2 J g2(p _2A;gcos^ + g2)2(p _2fcgcos^ + 2g2) ^ ' 

(1 - cos^ 9) [k^{d-l)- 2kqd cos 9 + 2q'^{d + 2 cos^ 9 - 2)] q'^'^dq sin'^-^ 9d9. 
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The integration over the angle variable cos^ can be performed explicitly. With the calcu- 
lations presented in Appendix, we get the effective pair correlator for n = 2, = 3 in the 

large-scale limit = | — > j : 

ctjr-\k ^ 0) = l^W0^!4^^2. (39) 
^ ' 40 ^ ' 

4 Energy dissipation and energy transfer 

The energy dissipation rate per unit of mass of an incompressible viscous fluid is given by 
the Navier-Stokes equations 



e = -zy y u(a;)Au(x)d'^x = p j ci'^x(Vu) 
Using wavelet decomposition ® for velocity field we get 



where 



e = -u I i2(ai, 02, bi - bajMai (01)^02(62) , (40) 



A)7k^^(a,k)^(-a2k)e-(--)|^ 



is the dissipation connection for the scale components. 

In symbolic form, the contribution of the fluctuations of all scales to the mean dissi- 
pation of energy per unit of mass can be written as 

e = ^z/ij / Ua,{xi)ua^{xj)d'^y.i(fyij, (41) 

where is the viscosity between and scales. Being well localized in both the real and 
the wavenumber space, the analyzing function perceives the interaction of the components 
of the same or close scales stronger than the contributions of the significantly different ones 
I log(ai/a2)| ^ 1. For the Daubechies wavelets (the orthogonal wavelets with compact sup- 
port) often used for numerical simulation of turbulence, the viscosity connection coefficients 
can be found elsewhere |34j . 
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For a qualitative estimation of the behavior of the viscosity connection as a function 
of the scale ratio ai/a2 of the interacting scale components, let us consider the vanishing 
momenta wavelet family of gaussian wavelets ()36|) . considered by Lewalle in a wavelet-based 
analysis of energy dissipation |T31E]. The viscosity connection ()4ip can be then evaluated 
analytically: 



J_^' ' (27r) 



-ir2 



(42) 



The main contribution to energy dissipation comes from the terms with coinciding or 
closed arguments x = 61 — 62 ~ 0. In this limit (with d = 1, taken for simplicity) we get 



2n + l 
2 



^n~^^f^^) (2n + l)!!. 

Introducing the ratio t = ^, we can study the behavior of the viscosity connection as a 
function of scale ratio 

2n + l 

The plot of the viscosity connection fi„(ai,a2,0) as a function of the scale ratio t = ^, 
for the first three wavelets {n = 1,2,3) of the (jHUI) family, is shown in Fig. El As it can be 
seen, regardless the number of vanishing momenta n, the dissipation term has a maximum at 
coinciding scales {t ^ 1). For this reason, if a discrete wavelet decomposition is used instead 
of continuous one, it is sufficient to keep the two main terms in the energy dissipation: the 
equal-scale interaction and the neighboring scale interaction 

e ului^ J dxijl{x)A^i^{x) - 2ul'\i^ j dxipf\x)Aipi^{x) + .... (43) 

The first term is the standard viscosity term, the last is the Kraichnan nearest scales inter- 
action. 

The nonlinear energy transfer between neighboring scales can be evaluated by consid- 
ering the wavelet connections corresponding to nonlinear term (uV)m of the Navier-Stokes 
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0.5 1 



Figure 3: Graph of the function ^i(^^_^_^2^ (f+t) ^ for = 1,2,3 related to the viscosity 
connection for gaussian type wavelets. 

equation. To keep with the wavelet turbulence cascade models jl], we restrict ourselves with 
a discrete wavelet transform with binary scale step, ao = ^: 



u(x) = u^Vkl^) + Error term, '^/'kl^) 



boal 



(44) 



Without loss of generality the mesh size is set to unity bo = l. In our consideration of the 
DWT representation ()44|) applied to hydrodynamic turbulence, in contrast to many schemes 
applied for numerical simulation of the NSE [HHl ITU] , we have no a priory arguments to 
assume a mutual orthogonahty of the basis functions ipl- To keep the wavelet decomposition 
fl44|) unique, the orthogonality of basic functions is not required. It is sufficient if the set of 
basic functions forms a frame, i.e. V/ G L^(]R), 3 A > 0, B < oo such that 



< B 



jk 



If A = B, the frame is called a tight frame. 

Assuming the basic functions ipl form a frame, and restricting ourselves with the case of 
incompressible fluid, we cast the NSE system in the form 



dxf^ 



1 dp'f^ipiix) 
p dx^ 



+ vu^^Aipiix), 



(45) 
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with Greek letters used for the coordinate indices, the bold face for the vector subscripts 
being dropped, and the summation over all pair-matching indices assumed. The component 
fields ul°' = are the functions of time only; so we deal with a typical cascade model. 

Our goal at this point is to derive the energy transfer between the components of j-th 
scale and the next small (j + l)-th one. Let us define the energy of the j-th scale pulsations 
as 

= ^ E Km = I d'xi,i{xWm{xl (46) 

ak 

with the assumed unit normalization of the basic function J d'^x\ip{x)\^ = 1. The contribution 
of the nonlinear term of the NSE to the time derivative of Ej is 

AE, = -AtuTu^iur J d'x,lf:,{xWm{xf-^. (47) 

For the orthogonal wavelets, that are most often used in numerical simulations only the 
terms of coinciding scales r = I = j survive in the r.h.s. of (@7j). The energy flux from the 
j-th scale to the next (j+l)-th scale is then proportional to \u^^/{boal), in exact accordance 
to the Kolmogorov phenomenological theory [Tj. In more general terms of nonorthogonal 
basic functions, the next term in the r.h.s. of (jTTj) is proportional to u^'^^u^^^uK This term 
can be interpreted as the material derivative u^'S/iu^'^'^Y of the mean energy ^"^^ ^ of the 
small scale fluctuations travelling along the stream of large scale velocity uK 

The energy transfer terms analogues to ()47|) have been already considered in the or- 
thogonal wavelet formalism by CMeneveau 0. They can be obtained directly from the 
component equations, by multiplying (PT|) by Uaiik). This leads to the energy transfer in 
(a, k) space 

j .i»(A:)il^l-^(k,q,k-q)<,(g)«,,,(fc-g)-^-^^^-^ (48) 

given in jBj. 

5 Kolmogorov hypotheses 

The Kolmogorov theory of the locally isotropic turbulence is formulated in terms of relative 
velocities 

5u{r,l) =u{r + 1) -u{r). (49) 

The probability distribution of relative velocities (jlUj) hardly can be studied by the Fourier 
transform, in case the velocity field u{r) is not homogeneous. According to Kolmogorov 
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PP, the turbulence in space-time domain G, is referred to as a stationary turbulence if, 
for any fixed u{r, t), the distribution of relative velocities 6u{r, t) is stationary and isotropic. 
Physical assumptions on the locally isotropic turbulence were formulated in terms of the first 
and second Kolmogorov hypotheses, that reside on the definition of the Reynolds number. 
This is not a quite rigorous mathematical definition. We shall show that the Kolmogorov 
hypotheses are the statements about the behavior of the scale components of velocity field. 

First, we have to note that the definition of the Reynolds number is consistent within 
the multiscale framework. In fact, by definition 

Rei = ^, (50) 

where ui are said to be "pulsations of the scale the rigorous definition of those can be 
given by virtue of wavelet components considering ip as an apparatus function used to 
measure the pulsations. Going further, we find the wavelet components (0) to be identical 
to velocity increments (jlHI), in case ip is the Haar wavelet: 

{1 < X < 1/2 
-1 l/2<x<l. (51) 
otherwise 

Now let us consider the Kolmogorov hypotheses pQ: 

HI: The first hypothesis of similarity For the locally isotropic turbulence with high 
enough Re the PDFs for the relative velocities are uniquely determined by the viscosity 
u and the mean energy dissipation rate e. 

H2: The second hypothesis of similarity Under the same assumption as for HI the 

3 1 

turbulent flow is self-similar in small (but still / ^ ijie"i) scales in the sense that 

5u{r, XI) X^6u{r, I), Xe M+. (52) 

Using the Taylor frozen flow hypothesis, and therefore considering one-dimensional pulsations 
u{x), according to the definition of wavelet coefficients (jSj), with ip = h we get 

Ui{r) = I ( — j — j u{x)dx = I dth{t)u{lt + r) 



— oo 







1/2 ^1 

u{lt + r)dt- I u{lt + r)dt. 

Jl/2 
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For small values of / we can approximate 



1 113 
Mr) ~ ^u{r + -I) - -u{r + -I), 

or, taking into account the statistical homogeneity of the flow, we get 



uiir) ^^u{r)-^u{r + ^). (53) 
So, the power-law behavior (fH^ . viz 

Mr) = -^5M(r,^) 

is just a particular case of a local regularity of wavelet coefficients, with the Haar function 
(j^Tj) being used as a basic wavelet. As it was shown in general settings [32], the wavelet 
coefficients W^{a,x)[f] of a square-integrable function /(x), which has the Lipshitz-Holder 
exponent h at the point x = Xo, behave as |W^(a, ~ a'^ inside the cone |a; — Xo| < const 
for any admissible wavelet -0 which satisfies the regularity condition 

dx{l + \x\)\iIj{x)\ < oo. (54) 

The condition ()54|) is rather loose, and in physical settings one can always assume that it 
holds for any analyzing function used to measure the pulsations of scale /. So, the second 
Kolmogorov hypothesis can be formulated as follows: 

H2: Generalized second Kolmogorov hypothesis of similarity Under the same 

3 1 

assumption as for HI the turbulent flow is self-similar in small (but still / ^ u^e'i) scales 
in the sense, that the pulsations of the turbulent velocity deflned as 

ui{b) = y ("~r~) 

where ip{x) is any analyzing function satisfying the admissibility condition (jS)) and the reg- 
ularity condition (jSH), have the following power- law behavior 

M&)r = r h=l, (55) 

for all spatial points b occupied by turbulent media. 
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6 Concluding remarks 



In this paper we have a stochastic hydrodynamics approach to the NSE based on a wavelet 
decomposition of both the velocity field and the stirring force used to balance the energy 
dissipation. In spite of a good deal of papers devoted do different choices of the stirring 
force for the NSE in wavenumber space - see e.g. . IH7| ITT] and references therein - 
what is essentially new in our approach is the definition of random force by the correlation 
function of its wavelet coefficients. Establishing the force correlator in the space of wavelet 
coefficients, we have got it easy to get rid of loop divergences in the stochastic perturbation 
expansion, having at the same time desired physical properties of the forcing (energy injec- 
tion at a given scale). Therefore, a new UV-finite framework is constructed for statistical 
hydrodynamics. This is a technical framework for analytical evaluation of the statistical 
characteristics of turbulent fiuctuations, such as their correlation and response functions, 
by means of continuous wavelet transform. Deriving physical consequences, such as energy 
cascade between scales or scale-dependent corrections to response function, we specially did 
not touch the renormalization group aspects |15^ IHS] of the problem and multifractal formal- 
ism jSnmniEl- in fact, both are related. The former is a generalization of the description 
of hydrodynamic turbulence in terms of differential equations to the description in measure 
settings. Partially, the relation of wavelets and RG in turbulence description is discussed in 
|41j . and will be studied in more details in connection to multifractal properties of hydro- 
dynamic turbulence. Besides, the possible comparison with the RG based classification of 
asymptotic regimes of isotropic turbulence f7| |^ can be considered as another perspective 
of the proposed method. 

The study of hydrodynamic turbulence by methods of stochastic differential equations 
and those of quantum field theory has at least half a century history. Regardless phenomeno- 
logically clear and widely accepted Kolmogorov P theory of fully developed turbulence, still 
there are discussions on the preference of either differential equations, or field theoretic meth- 
ods based on renormalization group, or multifractal approach to describe the turbulence in 
an incompressible fiuid. 

As it concerns the physical interpretation of the velocity field wavelet coefficients, by 
this paper we intended to say that stochastic nature of spatially extended hydrodynamic 
turbulence prescribes a certain kind of "wave-particle dualism" to the turbulent phenomena, 
in a sense, that the answer we get depends on the basic functions used to describe the 
turbulence. If the basis of plain waves was chosen, there are no fair reasons to comply about 
k^O behavior or paradoxes: what we get is what we set. Alternatively, if we want to have an 
analytical description of the spatially extended turbulence compatible with the Kolmogorov 
phenomenology of local turbulent pulsations of given scales, we need to set a functional basis 
that respects the scale locally. This is the wavelet decomposition. 
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The Fourier transform, being essentially nonlocal, apparently does not fit the above men- 
tioned requirements, but the windowed Fourier transform, or wave packet decomposition used 
by V.Zimin |3] and T.Nakano for the analysis of turbulence does, and possibly there is only a 
technical difference between our approach and that of Nakano [5^ . However, it is important 
to emphasize that the incorporation of the basic function ip into consideration makes us to 
admit that the definition of the local fluctuations of a given scale is not completely objective 
and depends on means of observation. As it was shown, the Kolmogorov hypotheses (K41) 
were easily rewritten in the wavelet framework for a multiscale description of turbulence. 
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A Calculation of one-loop diagrams 
A.l Response function 

To calculate the one-loop diagram in the response function, we introduce the following tensor 
structure 

L{k, q, a, s) = —m{k, a, b, c)o{q, b, l)m{p^ c, /, s), (56) 

where the summation over all dummy indices is assumed. The following notation for the 
orthogonal projector and the vertex (@]) is used {p = k — q): 

o{q, b, I) = 6bi - m{p, c, /, s) = ^ [pio{p, c, s) + Pso{p, c, /)] . 

After all convolutions in matching pairs of indices, substituting k ■ q = kqfi, where /x = cos 9 
is the cosine of the angle between k and q, we get 



4 ■ — 4p2 
-A;2/i^ + A;g/i(2/i2 - 1) + /iV , -2k^fx + 2k'^q{2fx'^ - 1) + Skfxp'^ 

2qp'^ 
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After substitution = k'^ — 2kqfi + in the numerators and some algebraic simplification 



L{k, q, a, s) = ^ [Tik^6as + T2kaks + T^^kaQs + T^qaks + T^qaQs] 



where 



T2 



cos^ e - 1 

2 ' 

A;^ + g2 - 2 cos^ 9{k'^ + 2g2) + Akq cos^ 6' 



2p2 

cos^ 6 — kq cos 6^ 



3 9 ) 

k^ cos 6* — 2k'^q cos^ 6* + S/cg^ cos — 2k'^q 
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2qp^ 

k'^ — kq cos 9 



5 



p2 



To calculate the whole one-loop integral contribution to the response function, the tensor 
structure L{k, q, a, s) is multiplied by the integral over the frequency component 

£ '^MG.i. - ,) ^ (57) 

Using this structure we easily get the IPI one-loop contribution to the response function 



u^(k) = Go{k)ri^{k) + Gl{k)AX^ j q'^-^dqde sin'^-^ 



A(q) 



2z/2q2 fco + q2 + p2 



^ [Tik^6as + T2kaks + T^kaqs + T^qaks + T.^qaqs]fls{.k). 

In (i = 3 we always assume k = fce^, with 6 being the polar angle k ■ q = kq cos 6 and (j) 
being the azimuthal angle. Let us evaluate the integral (jSHj) in ci = 3 

u,ik) = Go{k)Uk) + \^Gl{k)S2 I q^dqdOsme ] ^, ^^f ^ 

J v^fT ^ + + 

diag(Tifc2 + ^ sin^ ^Ts, Tifc^ + ^ sin^ ^Tg, 0)] Uk) 

Here diag() means diagonal matrix, where the first two terms give the transversal contribu- 
tion to viscosity, and the last term, giving the longtitudal contribution, vanishes identically: 

fc2(Ti + T2) + kq cos ^(Tg + T4) + q^ cos^ dT^ = 0. 
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Finally, after integration in angle variable cos 9, we get 

o roo 

i{k) = Goimik) + X'Glik)^ / dqA{q)dmg{R{k/q), R{k/q), Q)Uk) (58) 

^ Jo 



Ui 

where 

1 r 

12x - 16x^ - 8x^ + (x^ - 1)^ In 



R{x) ^ 



32x3 



1 \ 2 
X — 1 



X+1 



+ (3x^ - 2x^ + 12x2 -8) In ^^l^^^l. (59) 
A. 2 Correlation function 

Similarly, for the tensor structure of the one-loop contribution to the response function, we 
evaluate the tensor structure of the one-loop contribution to the pair correlator 

C{k, q, a, s) — m{k, a, 6, c)m{k, s, I, t)o{q, b, l)o{p, c, t). 

After algebraic simplification 

-k"^ + 2k^fiq + fc2(-/ip2 - ii'^q^ + 2p^) 



C{k,q,a,s) 



-\- kak. 



4p2 

k\l + fi^) - 2kqii{l + 2/^2) + 2(2/x^g2 _ ^2) 



+ {kaqs + qaks) 



4p2 

-k'^li + Ak'^qyi^ + kii{-Aii'^q^ + q^) 



Ap^q 

- Ak^iiq + k\A^^q^ - - q^) 
Ap^q^ 



ft, - ^A, -r ft. y^-tp^ q - y — q /^^x 

+ ^aQs ■ (60) 



The trace of this tensor structure, i.e. C2{k,q) — '^^C{k,q,a,a), required for the energy 
spectra evaluation, is equal to 

C2(A;, q) = [k^d - 1) - 2kqpid + 2q\d + 2//^ - 2)] . (61) 

In the important case of c? = 3 this gives 
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The integral in frequency argument of the product of squared response functions in the 
integral (j34j) gives in the limit of zero frequency /cq— ^0: 



y_oo 27r 2z/'^ q^(k - q)^ + (k - q)^ 



(63) 



In the important case of d = 3, with the Mexican hat wavelet (n = 2) taken for defi- 
niteness, the integral over the angle variable = cos 9 can be evaluated analytically for the 
single band random force ()25p . The angle integration in (jHHj) gives 



C. 



d=3,n=2 
e// 



{k) = X- 



,k^\Go{k)\H7ralDl 
2^ 



8 n2 /-oo 



dqq\cik/q), 



(64) 



where 



2^ X 



2 + — 2x?/ 



F(a;) + F(-a;) 



(65) 



g-a2g2(2+x2-2x) 



(-3 - 2alq^{-l - 3x + x^) + 2a^q\-l - 2x - 2x2 + 2x^) 



+ 2a^g^(2 + 2x + x^)) + 2(4 + x^)ExpIntEi {-alq'^{2 + x^ - 2x)) 



(66) 



In the limit of small wave numbers x = > this gives 

q o 



dqq lim Ux) = - 



Therefore the effective correlator tends to 



^d=3,n=2 
^eff 



(fc 0) = 



7 A;2|Go(A;)p7rfa3D2 



40 



Z/' 
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